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Abstract 



We present results from a lattice hadron spectrum calculation using three fla- 
vors of dynamical quarks — two light and one strange, and quenched simula- 
tions for comparison. These simulations were done using a one-loop Symanzik 
improved gauge action and an improved Kogut-Susskind quark action. The 
lattice spacings, and hence also the physical volumes, were tuned to be the 
same in all the runs to better expose differences due to flavor number. Lattice 
spacings were tuned using the static quark potential, so as a byproduct we 
obtain updated results for the effect of sea quarks on the static quark po- 
tential. We find indications that the full QCD meson spectrum is in better 
agreement with experiment than the quenched spectrum. For the ++ (ao) 
meson we see a coupling to two pseudoscalar mesons, or a meson decay on 
the lattice. 
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I. INTRODUCTION 



Computation of the properties of hadrons beginning from the QCD Lagrangian is a 
major goal of lattice gauge theory, and steady progress has been made. The computational 
burden of including dynamical quarks is a major obstacle in the use of lattice QCD to 
compute hadronic properties. As a result, many quantities are much better determined in 
the quenched approximation than in full QCD. One way of studying the effects of dynamical 
quarks is to calculate in quenched and full QCD, using the same valence quark action in 
both cases, and matching the lattice spacings and physical sizes of the lattices, so that any 
differences that are found can convincingly be ascribed to the dynamical quarks. Here we 
present a calculation of the hadron spectrum at a lattice spacing of about a = 0.13 fm, using 
quenched and full QCD lattices at the same lattice spacing. The lattice spacing was tuned 
by making short runs on smaller lattices, adjusting the parameters to match the lattice 
spacing of an initial quenched run at 10 /g 2 = 8.0. 

We use an improved Kogut-Susskind quark action which removes tree level order a 2 
lattice artifacts 0. The gauge action is a one-loop Symanzik improved action ||. This action 
has been shown to reduce flavor symmetry breaking and to improve rotational symmetry of 
the hadron spectrum, and to give improved scaling of hadron masses as a function of lattice 
spacing ||. 

Another important improvement of these calculations over previous generations is that 
we use three flavors of dynamical quarks. For quark masses larger than the strange quark 
mass we use three degenerate flavors, and for light quark masses less than m s we use two light 
flavors, keeping the third quark mass at about the strange quark mass. We have also done 
one two-flavor simulation on a matched lattice to check for effects of the dynamical strange 
quark. For the runs with 2 + 1 dynamical flavors, we computed hadron spectra using valence 
quark masses equal to the sea quark masses. In the quenched run we computed hadron 
masses with valence quarks with the same masses, and nondegenerate propagators using a 
strange quark mass of am s = 0.05. Finally, in the two dynamical flavor run we computed 
hadron propagators using light quark masses equal to the dynamical mass, am u ^ = 0.02 
and a non-dynamical strange quark with am s = 0.05. 

Two issues that we do not completely address are extrapolation to zero lattice spacing 
and extrapolation to the physical light quark mass. Using matched quenched and full QCD 
lattices allows us to draw conclusions about the effects of dynamical quarks without explicit 
continuum extrapolation. Although it is in principle possible that the discretization errors 
in quenched and full QCD are very different, we expect that the differences in lattice spacing 
dependence are in fact comparable to the differences in physical quantities themselves. Since 
corrections to scaling are in any case quite small with our current improved action pi|3[] (see 
also below), we confidently expect conclusions drawn at fixed lattice spacing to survive in 
the continuum limit. We are beginning a series of simulations at a smaller lattice spacing 
which will eventually allow us to make the continuum extrapolation. A few preliminary 
quenched points from these finer lattice runs are included here to give an idea of the size 
of these effects. A complete chiral extrapolation will be more difficult. In this work, we 
attempt an explicit chiral extrapolation only for the shape of the static quark potential, and 
show other quantities as functions of the quark mass. Some quantities such as "J" f|] are 
only minimally sensitive to chiral extrapolation, and such quantites provide immediately 
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useful tests of dynamical quark effects. 

In addition to presenting hadron spectra, we update our computation of the static quark 
potential. For this quantity we have sufficiently accurate data that we can hazard an extrap- 
olation to zero quark mass to produce numbers that can be compared with phenomenological 
potential models. Because the static potential is determined very accurately, it clearly shows 
differences between quenched and full QCD. In fact, one can even see the differences be- 
tween two and three dynamical flavors, and a "kink" in the plots where we change from 
three degenerate flavors to two light and one heavy indicates a noticeable difference between 
two light plus one heavy and three light dynamical flavors. 

In the meson sector we find differences between full and quenched QCD. A nice way 
of exposing these differences is the ratio J proposed by Lacock and Michael 0. We find 
that this quantity is in better agreement with experiment in full QCD than in the quenched 
approximation, as predicted in ||. This is consistent with results of the CP-PACS || and 
JLQCD [[]] collaborations, who also conclude that inclusion of two flavors of dynamical 
quarks improves agreement of the lattice spectrum with the real world. 

In the isovector J p = ++ (ao) channel we find a large difference between quenched 
and three-flavor results. We interpret the three flavor results as an avoided level crossing 
between the ++ meson and a two pseudoscalar state. In other words, we see the a decay 
to two mesons. 

We include tables of the mass fits we have chosen, so the reader can compute his or her 
own favorite mass ratios. 

II. SIMULATION PARAMETERS 

For our two and three flavor simulations we used the standard hybrid-molecular dynamics 
"R algorithm", |7j with one pseudofermion field for runs with degenerate quarks, and two 
pseudofermion fields for runs with different strange and up and down quark masses. In all 
cases, we used trajectories with unit length in the simulation time. Basic parameters of 
these runs are summarized in Table |I[ 

Two sources of systematic error in this method are the accuracy of the approximate 
sparse matrix solution required at each time step in the integration of the molecular dy- 
namics equations, and the effect of the nonzero time step used in the integration. We have 
investigated these effects in the 2+1 flavor simulations at 10/ g 2 = 6.8 with two flavors of 
quarks at mass am u ^ = 0.02 and one flavor with am s = 0.05 on a 16 3 x 48 lattice. (These 
masses are approximately 0.4 times the physical strange quark mass and the physical strange 
quark mass, respectively). Figures |Tj and ^| show the plaquette and ipip as a function of the 
conjugate gradient residual used in the updating. Finally, to see how this effect propagates 
into the hadron masses, we show the Goldstone pion mass in these same runs in Fig. |[ Since 
the effect of this error is poorly understood, we adopted a conservative choice of 1 x 10~ 4 or 
5 x 10 -5 in most of our runs, and used 2 x 10~ 5 for the heaviest quark mass (am q = 0.40) 
where convergence was very fast. 

The effect of integration step size is better understood. We verified the expected 
quadratic dependence of the error on the step size in our three flavor code. The first panel 
in Fig. |] shows the result of one such test, where we ran on a quite coarse lattice using the 



4 



conventional action, comparing the old code with three degenerate flavors to the two-plus- 
one flavor code with m u ^ = m s = 0.02/a. (The "old code" uses a single pseudofermion 
vector, with a weight of 3/4 in the fermion force, while the "two-plus-one flavor" code uses 
two pseudofermion vectors, one with a weight of 2/4 and the other with weight 1/4. In this 
test, both pseudofermion vectors had the same mass. In each case, the multiplication of 
the Gaussian random vector by M* to produce the pseudofermion vector was done at the 
appropriate point in the time step to make the error in a single time step order e 3 , making 
the accumulated error over a trajectory order e 2 . It can be seen that both sets of points 
extrapolate to the same limit, although the size of the effect is quite different. The right 
hand panel of Fig. |] shows a similar plot from one of our pre-production tuning runs with 
the improved gauge and quark action, at 10 /g 2 = 6.80 and am q = 0.05, on a 16 3 x 48 lattice. 
This is close to the value of 10/ g 2 used in the production run at this quark mass. Note that 
the finite step size corrections are quite sensitive to the action being used - even the sign 
of the effect differs between the two actions. Based on these tests and previous experience, 
we used a step size of no more than two-thirds of the light quark mass or 0.03, whichever 
was smaller, in our production runs. In the production run at am q = 0.05 we used a step 
size of e = 0.02. From the slope in Fig. |] we can infer that this caused a fractional error 
in (tjjip) of about 0.004, and a similar analysis for the plaquette suggests a fractional error 
of about 0.0006. We also looked at pion masses and some of the Wilson loops involved in 
the computation of the static quark potential, and we were unable to resolve statistically 
significant effects on these quantities in our tuning runs. 

In our production runs with dynamical quarks we measured the potential and the spec- 
trum at intervals of six simulation time units, and archived these lattices. The autocorre- 
lation of the plaquette at six time units, or successive measurements, was generally about 
0.1. We investigated the effect of autocorrelations on the potential and spectrum by block- 
ing together different numbers of measurements before doing the fitting. For the potential 
measurements we chose to block five measurments (30 time units). For the hadron spectrum 
most of the particles showed no systematic effects of blocking measurements together. The 
exception was the pseudoscalar mesons, where blocking gave significantly better confidence 
levels and larger error estimates. For the pseudoscalars we chose to block together four 
measurements, or 24 time units. 

The physical size of our lattices is La = 20 a « 2.6 fm, which is similar to or larger 
than other recent full QCD simulations. Basic parameters of several of these calculations 
are summarized in Ref. ||. 
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FIG. 1. The effect of the conjugate gradient error used in the updating on the plaquette in a 
three flavor run with quark masses 0Am s and m s . 
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FIG. 2. The effect of the conjugate gradient error used in the updating on ipip in the same 
three flavor run. 
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FIG. 3. The effect of the conjugate gradient error used in the updating on the Goldstone pion 
mass in the same three flavor run. 
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FIG. 4. (i/iip) versus the squared step size. The left hand panel is an algorithm test done on 
a 12 lattice using the one-plaquette gauge action and conventional quark action at 6/g 2 = 5.10 
with three quark flavors with mass am q = 0.02. The octagons use one pseudofermion field with a 
factor of 3/4 in the force term, appropriate for three flavors, while the squares use the 2 + 1 flavor 
code, with separate fermion force terms for one and two flavors, but with the same mass for both 
terms. The right hand panel shows {^tp) using improved gauge and quark actions at a ~ 0.14 fm. 
(10/g 2 = 6.80 and am Ui d = am s = 0.05). 

III. LENGTH SCALES FROM THE STATIC POTENTIAL 

The static quark potential is often used to define the length scale in lattice simulations. 
Advantages of using the potential include ease and accuracy of its computation, and its 
lack of dependence on the valence quark mass. In comparing quenched and full simulations, 
subtleties arise because the potential does depend on the masses of the sea quarks. In 
Ref. |§ we demonstrated the effects of sea quarks on the potential using this improved action. 
Because these effects are important in our analysis of the hadron spectrum, we update and 
extend these results here. Our methods for computing the potential and our reasons for 
using rq, a variant of the conventional r ||10|| , are described in ||. r is conventionally 
denned by r^F(ro) = 1.65, and rq by rfF(ri) = 1.00. 

The fitting form used here is slightly more complicated than the form used in Ref. ||, 
with an extra term to take into account lattice artifacts at the shortest distances. Following 
a procedure used in Ref. JTTj ] , 

V(r) = C + ar-a/r + X (V free (f) - 1/r) . (1) 

The last term, used for r < 2.5, approximately compensates for remaining lattice artifacts. 



S 



Here Vf ree (r) is the potential calculated in free field theory, using the improved gauge action. 
Adding this term to the fits significantly improves the goodness of fit and makes the fit 
parameters less sensitive to the choice of distance range. For the a »s 0.13 fm runs we 
typically find A « 0.3 - 0.4. 

In Figs. H and || we show the dimensionless quantities r y/a and r\y/a respectively as 
functions of the quark mass, represented by {m^/nip) 2 . This places the quenched approxi- 
mation at (m^/rrip) 2 = 1, and the chiral limit at the left side of the graph. In these plots 
the octagons are runs with three degenerate sea quarks, except for the rightmost point 
which is the quenched limit. Squares are runs with am s = 0.05, its approximate physical 
value, and am u ^ < 0.05. The isolated diamond is our two flavor run. Finally, the cross 
at (m^/mp) 2 = 1 is the finer lattice quenched run. From the two quenched points we see 
that remaining lattice artifacts are small compared with the effects of the sea quarks. In 
particular, the central values for r^^fo and ry^fa changed by less than 1% when the lattice 
spacing was reduced from 0.13 fm to 0.09 fm, a change of about 35%. The kink in the plots 
at (tott/top) 2 ~ 0.46 (am u> d = 0.05) shows the transition between three degenerate flavors 
and "2 + 1" flavors. We can clearly see the distinction between two and three flavors, as 
well as the effect of using two light and one heavy flavor rather than three degenerate flavors 
(the "kink' at (m w /m p ) 2 « 0.46). 

If we extrapolate r y/a to the physical quark mass, as shown in Fig. |, we find r Q y/a = 
1.114(4) (statistical errors only), a number which can be compared with phenomenological 
potential models. The two quenched points give some idea of the possible systematic error. 
Since the squared lattice spacing in the finer lattice is about one half that of the coarser, we 
might expect a shift of about twice the separation of these points in the continuum limit. 
Since the error bars on these points overlap, we don't know this systematic error well enough 
to justify such an extrapolation at this point. We expect that the effect of this systematic 
error is mostly an overall shift of the graphs, but the next generation of simulations should 
clarify this. 

While using r\ to define the length scale has the advantage that it can be done more 
accurately, tq has the advantage that it has been related to phenomenological potential 
models, which consistently place it around 0.5 fm |H|. Therefore, to estimate r% in physical 
units, we plot ro/ri in Fig. |f[ Extrapolating linearly in (m n /m p ) 2 to the physical value gives 
r o/ r i = 1-449(5) (statistical error only), or with tq = 0.5 fm, r\ = 0.35 fm. 
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FIG. 5. Effects of dynamical quarks on the shape of the potential. Here we plot r^^fo as a 
function of the quark mass. The two quenched points are at the far right, with the octagon coming 
from the 10 /g 2 = 8.0 run and the cross from the 10/ g = 8.4 run, which has a lattice spacing of 
about 0.09 fm. The remaining octagons are full QCD runs with three degenerate flavors, and the 
squares are full QCD runs with two light flavors and one heavy. The diamond is the two flavor 
run, and the burst at the left is a linear extrapolation of the 2 + 1 results to the physical value of 
(ra^/rap) 2 . 
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FIG. 6. The same as Fig. ||, except we plot r\y/a. Physically, the difference is that this 
quantity is sensitive to shorter distances than rQyfa. 
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IV. PROCEDURES FOR DETERMINING MASSES 



All of our hadron propagators used wall sources and local sink operators. Several different 
wall sources were used. For the "pointlike" hadrons, for which all the quarks can be on a 
single corner of the hypercube, a "corner wall" source gave the best results. This source 
is simply a 1 on each (0, 0, 0) corner of the 2 3 cubes on the chosen time slice. However, 
to isolate the decuplet baryons a non-pointlike source is essential, and we used "even" and 
"odd" wall sources, where 1 or (— iy + y +z is placed on each site, respectively. This set of 
sources was developed in Ref. ||12|| . For the nonlocal pseudoscalar and vector mesons we 
used two wall sources made from empirically determined linear combinations of the nonlocal 
pion operators. Finally, for the nonzero momentum mesons we used a quark source with 1 
on each site, and an antiquark source with e tk ' x on each site. All of the configurations were 
gauge fixed to the Coulomb gauge before computing the propagators. 

In most cases we computed propagators from four source times evenly spread through 
the lattice (only one source slice was "turned on" at a time). For the corner source we used 
eight source time slices for the light quark particles with am q < 0.04, and we used eight 
source time slices for the "even" and "odd" source baryon propagators. (The A propagator 
is very noisy, and propagators computed from source times separated by 8 x 0.13 fm were 
basically independent.) 

For Kogut-Susskind quarks the meson propagators have the generic form 

H(t) = A i ( e " mii + e"" 1 ^-*)) + ( e ~ Ki + e- m '^ Nt -*A . (2) 



Here the oscillating terms correspond to particles with opposite parity from the ordinary 
exponential terms. Baryon propagators are similar, but have antiperiodic boundary con- 
ditions and the "backwards" terms include an extra factor of (—1)'. In most cases only 
one mass with each parity is included in the fits, but for half of the pseudoscalar meson 
operators the opposite parity terms are not present, and for the P-wave mesons we found it 
necessary to keep two simple exponentials. The quantum numbers for the various operators 
are tabulated in Ref. [|i~3H . 

Hadron masses were determined from fits to propagators, using the full covariance matrix 
to estimate errors. The maximum time distance used in the fits was chosen to include points 
with fractional error less than 0.3. Because of the oscillating components in the staggered 
quark propagators, it sometimes happens that the fractional error exceeds the threshold at 
one distance but is smaller at larger distances, so the complete criterion for the maximum 
distance included is the largest distance such that the fractional error on each of the next 
two points exceeds 0.3. Since the points at largest distance contribute little information, 
the exact large distance cutoff is not critical. To choose the minimum distance included, we 
first went through the fits and chose a minimum distance for each hadron in each dataset, 
choosing a distance where the confidence level was reasonable and where the mass appeared 
to reach a plateau. As expected, some propagators had larger(smaller) fluctuations than 
other similar propagators, resulting in the choice of a larger (smaller) minimum distance. For 
the particles of greatest interest, to reduce this effect we then "smoothed" these minimum 
distances, requiring that the minimum distances be smooth functions of quark mass and 
be the same for the quenched and dynamical runs. The resulting minimum distances are 
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strongly dependent on the quark masses, with smaller minimum distances for smaller masses. 
Most of this is due to the larger statistical errors at smaller quark mass, which result in the 
excited state contributions disappearing into the noise at shorter distance. However, it is in 
part physical, since splittings between the ground and excited states are larger for smaller 
quark masses. Table [Tl] shows the minimum distances that we chose, and the number of 
particles with each parity. 



V. RESULTS 



A. Pseudoscalar mesons 

We calculated propagators for all eight flavor combinations of the staggered quark pseu- 
doscalar mesons. These masses obey the "partial flavor symmetry restoration" predicted by 
Lee and Sharpe [[L4[] to very good accuracy. Specifically, Ref. |14| predicts that the leading 



order flavor symmetry breaking effects, which are order a 2 , leave degeneracies between pairs 
of pseudoscalar mesons for which 7 is replaced by 7^ in the flavor structure. For example, 
the local non-Goldstone pion, 7o75<S>7o75 in the "spin ® flavor" notation, is degenerate with 
the distance one pion, 75 ® 7^75, to this order. Moreover, all of the squared pion masses 
should depend linearly on the quark mass with the same slope to lowest order. Figure |8| 
shows this behavior for the quenched pion masses for am u ^ < 0.05. The results for the full 
QCD runs are similar. Since the local 75 <8> 75 pseudoscalar has the correct chiral behavior 
(and the best scaling behavior), we will use this pseudoscalar in the rest of the analysis 
unless we specifically indicate otherwise. Results for the full QCD runs are similar, but the 
flavor symmetry breaking is somewhat larger. 

In Fig. H the relation between the squared pseudoscalar mass and the quark mass is 
clearly nearly linear. The deviations from linearity and the effect of the dynamical quarks 
can be exposed by plotting the squared pseudoscalar masses divided by the quark mass, 
in Fig. This is essentially (ijjijj)ri/ f% with an (unknown) renormalization factor. This 
plot contains pseudoscalar mesons with both light and strange valence quarks (pions, kaons 
and "unmixed ss's"). There is clearly a systematic difference between quenched and full 
QCD. This difference increases with decreasing quark mass, and the two-flavor point falls 
in between the quenched and three-flavor points. The bursts among the quenched points 
are from the 10/g 2 = 8.4, a ~ 0.09 fm run, showing gratifying agreement with the a ~ 0.13 
fm points. Unfortunately, a coarser three flavor lattice, a ~ 0.2 fm, shows a large effect, so 
we would not want to use much coarser lattices in studying this effect. We do note that 
we expect scaling violations to be similar for the quenched and dynamical theories, so it is 
an advantage to have runs with matched lattice spacings. The deviations from linearity of 
m 2 are similar in quenched and full QCD. The upturn for larger quark masses signals the 
beginning of the transition to the heavy quark regime, where m 2 ~ m 2 . We do not fully 
understand the shape of this plot for small quark mass. There are several ways to interpret 
the difference between quenched and dynamical results. One could say that (ipip) / / 2 is too 
small in the quenched approximation, or one could say that the quark mass at which a 
desired m^/rrip is reached is larger in the quenched approximation than in full QCD. This 
second interpretation is consistent with CP-PACS results on the quark masses, in which 
they find that the quark mass needed to reach a given value of m v /m p is smaller in two 
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flavor QCD than in quenched QCD ||15|| . (Indeed, one could even use this quantity as a 



length scale, and conclude that r x is different in quenched and full QCD.) 

The largest part of the error bars in Fig. [5] come from the uncertainty in r\. However, this 
uncertainty is common to all of the points coming from the same set of lattices. In particular, 
all the a ~ 0.13 fm quenched points are correlated in this respect, as are the three two-flavor 
points. If we are interested in the dependence of the pseudoscalar mass on the quark mass 
on a fixed lattice, we may want to consider only the error from the determination of the 
meson mass in units of a. The left hand panel in Fig. [10] shows m\r\j{r\{m\ + m 2 )) for the 
quenched calculation, including only the error from amps and showing only the reasonably 
light mass points. In this panel the octagons are "pions" , with mi = m 2 = am u ^ and the 
bursts are "kaons", with arri2 fixed at 0.05, which is approximately the physical value of the 
strange quark mass. We see that this quantity is dependent only on the sum of the quark 
masses to very good accuracy. The center panel contains the same plot for the three-flavor 
runs, where now r\ is determined independently in each run. In this panel the octagons are 
runs with am u ^ = am s > 0.05, the squares the run with am u ^ = 0.04, the diamonds from 
an^u,d = 0.03, the crosses from am U) d = 0.02 and the bursts from am u ^ = 0.01. The three 
symbols for each of the runs with am U) d < 0.05 correspond to the "pion" , with both valence 
quarks light, the "kaon", with one light and one strange valence quark, and an "unmixed 
ss", with two valence quarks of mass approximately equal to that of the strange quark, 
but no qq annihilation. This graph is far from smooth, but most of the scatter comes from 
the fact that each set of dynamical quark masses has an independent uncertainty in r%. 
Note that each dynamical run (for example, the three bursts) shows qualitatively the same 
behavior as the quenched case, with the light-light pseudoscalar tending to a larger value. 
An interesting question is how the sea quark mass affects the pseudoscalar mass. This can 
be investigated by looking at the points at a(rrii + m 2 ) =0.1, which are the "unmixed ss 



points with both valence quark masses equal to 0.05/a. In the right hand panel of Fig. [LG 
we plot these points as a function of the light quark mass. There is a noticeable effect, with 
smaller light quark mass producing larger ss mass. The direction of this effect is consistent 
with the smaller pseudoscalar masses in the two flavor and quenched calculations seen in 
Fig. |. 

The selected pseudoscalar meson mass fits in units of the lattice spacing are tabulated 
in table [TTI]. In addition to the "pions" , the table also contains fits with one quark at about 
the strange quark mass and one lighter quark, or "kaons" . For the two and three flavor runs 
we also tabulate "unmixed ss mesons", with two valence quarks with am v = 0.05. 

We have attempted to fit the results in Fig. |9]to the forms predicted by chiral perturbation 
theory. In the quenched case, the behavior of m\ as a function of quark mass is derived in 
Refs. [[Ulpj]]. We use eqn. (9) in fTB] , with the parameter a, which is believed to be small, 
set equal to zero, and the analytic correction term added to the chiral log. For pions, then 
one has 

^ = C(l-51og(^) + Km), (3) 

where mi = mi = m, C and K are constants, and the chiral scale A may be taken as the 
r] mass. As shown in Fig. §, the fit to the a = .13 fm quenched data is good. It gives 
5 = 0.061(3) (statistical error only), which is on the low side but compatible with the range 
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reported by CP-PACS |L8| and is in excellent agreement with the result of Bardeen et al. 
I|, 5 = 0.065(13). 



Unfortunately, our attempts to fit the 3-flavor pion data in Fig. |9] to the corresponding 
full QCD chiral form [2(J have been unsuccessful to date. We can fit the 5 lowest-mass pions 
with reasonable confidence level, but the coefficients of the analytic terms are unreasonably 
large, and the fit misses the next lightest pion by a wide margin. If we try to fit the 6 
lowest-mass pions, the fit has terrible confidence level. Finally, good fits can be obtained 
by introducing, as an additional free parameter, an overall coefficient in front of the chiral 
logs. However, the value of that coefficient in the fit is much smaller than its predicted EO 
value. We are continuing to study this puzzling situation. Our current running at smaller 
lattice spacing may provide additional insight here. 
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FIG. 8. Squared pseudoscalar meson masses versus quark masses. These results are from 
the quenched runs. The octagons are the local pions (75 (g> 75 and 7075 <8> 7o7s), the diamonds the 
distance one pions, the bursts the distance two pions, and the squares the distance three pions. The 



degeneracies predicted in Ref. 14] are clearly visible. The lines are not fits; they simply connect 
the points. The column of pluses is from the quenched a ~ 0.09 fm run, showing the expected 
improvement in flavor symmetry with decreasing lattice spacing. Note that the Goldstone pion 
changes very little when the lattice spacing is decreased — the non-Goldstone pions come down to 
join it. 
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FIG. 9. The squared pseudoscalar mass divided by the quark mass in units of r±. The octagons 
are the quenched a ~ 0.13 fm run, and the squares the three flavor a ~ 0.13 fm runs. The diamonds 
are the single two-flavor run. The bursts are from the quenched a ~ 0.09 fm run. The fit is to the 
octagons (pions only) with r\(m\ + 777-2) < 1.4, using the form in eq. (||). 
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FIG. 10. Details of (mpsri) 2 /{rl(m\ + 012)). The three panels are the quenched results, the 
dynamical results, and the three flavor points with two strange valence quarks. 

B. Vector mesons 

We calculated propagators for the two local vector mesons, 7^ ® 7, (VT) and 7 7i (g> 7o7i 
(PV), and two distance one vector mesons, 7, ® 1 and 707^ ® 70- Any flavor symmetry 
breaking among these mesons is smaller than the statistical errors, so we simply quote 
results for the local 7$ <£>7i, or "VT" mesons. Table [TV] contains these masses in units of the 
lattice spacing. 



In Figure |TT] we plot the vector meson mass in units of rq versus the squared pion/rho 
mass ratio. In this plot there is a clear difference between the quenched and dynamical 
masses, with the full QCD vector mesons lying lower. However, the size of this effect 



depends on the length standard chosen, as illustrated in Fig. [T2|, where the same quantity 
is plotted in units of the string tension. Of course, the difference between these two plots 
simply arises from the difference in T\yfa plotted in Fig. 13. 
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a two-flavor result, and the squares three-flavor. The bursts are quenched a ~ 0.09 fm points, and 
the fancy plus an a ~ 0.2 fm three flavor point. 
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FIG. 12. The same as Fig. |lXJ, except in units of the string tension. 
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In our calculations we used a mass of am s = 0.05 for the strange quark mass. With the 
meson spectrum in hand, we can now go back and ask whether this choice was exactly correct. 
There are several quantities that we could choose to define the strange quark mass, and, 
especially in the quenched calculation, they will in general give different results. Perhaps 
the most accessible quantity is the ratio of the ss pseudoscalar to vector meson mass. (Our 
u ss pseudoscalar" does not include qq annihilation diagrams, so it is not the rj or rf.) We 



therefore tune this quantity using m sS = \J2m 2 K — m\ = 686 MeV and = 1020 MeV, or 
fnps/ m v = 0.673. 

In the quenched spectrum with am q = 0.05 we have mps/my = 0.643, indicating that 
am q = 0.05 is smaller than the desired strange quark mass. To compute the mass at which 
mps/mv = 0.673, we can do a linear fit to the vector meson mass as a function of quark 
mass, and combine that with a squared pseudoscalar mass proportional to the quark mass to 
conclude that the quenched strange quark mass defined by mps/my at this lattice spacing 
is am s = 0.058. In contrast, for the dynamical runs am s = 0.05 is a fairly good estimate of 
the strange quark mass, with m PS /m v taking values between 0.670(2) and 0.687(1). 

Differences between the quenched meson spectrum and the real world have been observed 
by the UKQCD collaboration |4[], and improvements of the spectrum when dynamical quarks 
are included have been reported by the CPPACS || and JLQCD || collaborations. In 
particular, the UKQCD collaboration studied the quantity 



dmv 

L PS 



J = mKt (4) 
ami 



where my and mps are the vector and pseudoscalar meson masses. This quantity has the 
advantage of being relatively insensitive to the quark masses, so that accurate tuning of the 
strange quark mass or extrapolation of the masses to the chiral limit is not essential. Of 
course, to compare to experiment the derivative in this expression must be replaced by a 
ratio of mass differences, and we choose 

(mg - m p ) 

J - mK *2{ m \- m l) ■ (5) 

Here m p is the mass of the vector meson including two light quarks, etc. We choose the 
m<£ — m p mass difference because the statistical error in J is dominated by the error in 
the vector meson masses, and the larger difference in — m p relative to, say, — m p 
leads to smaller statistical errors on J. Because all of the masses in Eq. [5] are correlated, 
we use a jackknife analysis to compute the error on J. Figure [13| shows the results for J 
in quenched and three flavor QCD. Following UKQCD, we plot this versus ms-./m^, for 
which the real world value is 1.8. The burst is the real world value of this definition of J 
(0.49), and the cross is the value of J found in the UKQCD quenched simulations. We see 
a clear effect of the sea quarks on this quantity. Indeed, any reasonable extrapolation of 
our data in rriK*/mK would pass near the real world point. Figure [13| also contains one 
point with two dynamical flavors. This point falls near the three flavor points, indicating 
that the dynamical strange quark is less important than the two light quarks. Although our 
quenched results are somewhat higher than the UKQCD value, they are significantly below 
the experimental value. The fact that the quenched points in this plot are to the right of 
the full QCD points is largely due to the fact that the mass of am q = 0.05 used for the 
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quenched strange quark was too small, as discussed above. If we use the observed slopes of 
the quenched vector meson masses and squared pseudoscalar masses as functions of quark 
mass to adjust these points to a quenched strange quark mass of am s = 0.058, the main 
effect is to shift the quenched points to the left. In particular, the rightmost quenched point 
moves to m/^/m^ = 1.78, but moves up by only 0.004. 
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FIG. 13. The mass ratio "J" in the quenched and full QCD calculations. Squares are the 
three flavor results, and octagons are the quenched results. The diamond is the two flavor run, 
using a non-dynamical strange quark with mass am q = 0.05. The burst is the real world value, 
and the cross is the UKQCD quenched value. The smaller error bar on the cross is the statistical 
error, and the larger the quoted systematic error. 



C. P-wave mesons 



Mesons with quantum numbers J PC = ++ , 1 ++ and l + ~, which are P-wave mesons 
in the nonrelativistic limit, are found as the oscillating "parity partners" of the 7r and p 
propagators with Kogut-Susskind quarks. Fitting these particles is difficult because of their 
larger mass. In order to get good fits, we need to allow two particles in the non-oscillating 
components, since we typically find that the excited 'V or "p" state is comparable in mass 
to the lowest oscillating (P-wave) state. Thus we find large errors, and must use small 
minimum distances to get good fits. In addition, plateaus in the effective mass are short. 
Figure [14| illustrates an example of a difficult but crucial case. This figure shows fits to 
the scalar (ao) mass as a function of the minimum distance included in the fit. With these 
caveats, selected fits for the P-wave mesons are shown in Fig. [L5] and Tables |V|, [VT] and 



VII, This figure and these tables contain several interesting features. For the oq and b\ the 



21 



full QCD runs give consistently better fits than the quenched run, and smaller masses for 
the lighter quarks. The b\ is consistently slightly heavier than the a%, although with the 
difficulties in extracting these particles we would not want to make too much of this. The 
diamond at the left in Fig. |15| is the experimental value for the a% and b\ masses. 

The scalar channel, "ao", is clearly very different in the quenched and full QCD runs. 
For large quark masses there is no visible difference, but as the quark mass is decreased the 
full QCD ++ mass drops below all the other masses. For all but the lowest quark mass, 
the quenched ++ is close to the other P-wave meson masses. We ascribe the behavior of 
the full QCD mass to the decay of the ao into tc + rj. (Bose symmetry plus isospin forbids 
decay into two pions.) Figure [16] illustrates this interpretation. In the figure we plot the 
quenched and full ++ masses versus quark mass. The straight line in the graph is a fit to 
the quenched mass for the heavier quarks, and represents the mass of a qq state. The curved 
line with the kink at am q = 0.05 represents the mass of 7r + 77. For am q > 0.05 we used three 
degenerate quark flavors, so the rj and ir are degenerate and this line is simply twice the 
pion mass. For am q < 0.05 we don't have direct information on the 77 mass, so we use the 
Gell-man-Okubo formula written in terms of an "unmixed ss" mass (just our pseudoscalar 
mass at am q = 0.05). 



m\= (ml + 2m 2 sS ) /3 (6) 
In the quenched case the ao mesons can couple to two-meson states through a "hairpin 



diagram" on one of the meson lines. Such diagrams, like Fig. 1(b) in Ref. [21], can behave 



like powers of t times e~ 2nint and therefore masquerade as a light ao when 2171^ < m ao . 
This may explain the lightest quark mass quenched point. In this analysis we used the a 
from the local source, or the "1 <g> 1" operator, which gave the best signal. Within the 
very large statistical errors, we saw relatively large breaking of flavor symmetry among the 
different ao channels, with some evidence that this reflects the masses of the different lattice 
pseudoscalars to which the various ao's should couple. 
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FIG. 14. Fitted masses for the scalar (ao) meson with quark mass am U; d = 0.03 as a function 
of the minimum distance included in the fit. In this plot the symbol size indicates the confidence 
level of the fit, with the symbol size used in the legend corresponding to 50%. Here the octagons 
are quenched results using two ^ (pion) states and one ++ (ao) state, while the bursts are 
quenched results with one state of each parity. The squares and diamonds are 2+1 flavor results 
using two and one h states respectively. 
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FIG. 15. P-wave meson masses in units of r±. The bursts are the quenched pseudovector 
mesons {a\ and b±), and the diamonds the full QCD pseudovectors. Where the difference is signif- 
icant the 1 H (pi) is heavier than the 1 ++ . The octagons are the quenched scalar (ao), and the 
squares the full QCD scalar. Crosses are the two flavor results. The diamond at the physical value 
of (mir/mp) 2 is the experimental a\ and b\ mass, and the two squares are experimental ++ (ao) 
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FIG. 16. ++ masses versus quark mass. The lighest fitted energy in the scalar channel. 
Octagons are quenched results, squares are three flavor results, and the burst is the two flavor run. 
The straight line is a crude extrapolation of the heavy quark points. The curved line is the ir + r/ 
mass estimate, as discussed in the text. The short vertical line marks the approximate quark mass 
where the ao mass is twice the quenched pion mass. 



We have evaluated propagators for baryons using the "corner wall" source for both 
degenerate and nondegenerate quarks, using a pointlike sink operator with all three quarks 
on the same lattice site. With nondegenerate quarks, the lightest states in this channel for 
zero, one or two strange quarks are the N, the A and the S respectively. However, since 
we took no special measures to make the operator with one strange quark and two light 
quarks orthogonal to the S, these propagators undoubtedly contain contamination from a 



D. Baryons 



25 



nearby E. In order to get the decuplet baryons, we followed Ref. |T2|], using a wall source 
on every spatial site, and using the operator in Eq. 6.3 of Ref. for the A. This is 
necessary because the corner wall source does not overlap this A operator. As a byproduct 
of the calculation of the decuplet mass, we obtain nucleon propagators from even site wall 
sources and a wall source containing all sites. These propagators are generally noisier than 
the corner wall source propagators, and the plateau in the effective mass occurs at larger 
distances. Therefore we generally use the corner source propagators. However, if one uses 
our fitting procedures on the even-wall or full-wall propagators, one invariably selects a 
smaller mass than from the corner-wall source propagators. This situation is illustrated by 
Fig. |T7|, which show mass fits to these two propagators as a function of the minimum distance 
included in the fit. In this figure the corner source propagators reach a plateau earlier and 
with smaller error bars. However, there are perfectly acceptable fits to the full wall source 
propagators giving masses significantly smaller than the corner source values. Perhaps the 
only good thing we can say about this situation is that the effect is similar in the quenched 
and dynamical runs, so as long as we are careful to make the same choices in both cases, we 
can investigate the effects of sea quarks on the spectrum while taking the statistical errors 
at face value. The results of the fits that we selected are listed in Tables |VIII| and [IX . 

The nucleon to rho mass ratio, or "Edinburgh plot", has long been used as a way of 
displaying lattice spectrum results. This ratio is known to be sensitive to lattice spacing, 
lattice volume and quark masses. As mentioned above, we do not address the issues of 
continuum extrapolation and chiral extrapolation in this paper. However, since our quenched 
and dynamical lattices are matched in lattice spacing and physical size, we are well positioned 
to ask if effects of dynamical quarks show up in this ratio. In Fig. ^ we show a variant of 
the Edinburgh plot, the "APE plot", where {m^/rrip) 2 is used as the abscissa, so that for 
small quark masses the abscissa is proportional to the quark mass. Most of the points on 
this plot are from the a ~ 0.13 fm matched lattice runs. It can be seen that there are no 
significant differences between the quenched and three flavor runs. The single two flavor 
point lies slightly above the trend, although this is probably not significant. This agreement 
between m^/m p for quenched and full QCD is in apparent conflict with our extrapolations 
of the conventional action ||23[| . The discrepancy, which may be due to residual discretization 
effects in one or both calculations, is under study. We are hopeful that currently running 
three flavor simulations at finer lattice spacing will shed new light here. Figure ^ also 
contains a point from a coarser lattice three flavor run and a preliminary point from a finer 
lattice quenched run (10/g 2 = 8.4, 28 3 x 96 lattice, a ~ 0.09 fm) These two points suggest 
that when we are in a position to do a continuum extrapolation the continuum results will 
be lower. 

It is interesting to compare these results to the conventional Kogut-Susskind quark ac- 
tion. Figure [19| shows the improved action results together with conventional action results 
at 6/g 2 = 5.7, 5.85. 6.15 [f24j] and 6.5 |2"5]| , which correspond to lattice spacings of about 
0.16, 0.12, 0.07 and 0.043 fm respectively. While a continuum extrapolation will be deferred 
until the 0.1 fm runs are completed, we can see in this plot that the improved action at 0.13 
fm gives results similar to the conventional action at a = 0.07 fm. 

Just as for the mesons, one of the biggest problems in comparing baryon masses to the 
real world is the need for an extrapolation in quark mass. For the mesons the quantity J has 
the nice feature that it is only minimally sensitive to this extrapolation. It is tempting to 
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try to construct similar quantities for the baryons. This suggests looking at the Q~ , the sss 
decuplet baryon. More generally, we could plot the mass of the decuplet baryon as a function 
of {mps/my) 2 to produce a variant of the Edinburgh plot which has the nice feature that 
there are two experimental points, one for ■m^-jm^ and another for m^/m p . Unfortunately, 
because the decuplet masses are hard to determine on the lattice, the errors are rather large. 
The result of this exercise is in Fig. |2(| Although the error bars and the scatter among the 
points are large, the overall trend of this plot is encouraging. It is premature to say whether 
this plot shows real differences between full and quenched QCD. 
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FIG. 17. Mass fits to nucleoli propagators as a function of the minimum distance included 
in the fit. These fits are from the three flavor run with 10/V = 6.85 and arriq = 0.05, and are 
fairly typical. Here the octagons are from propagators with a "corner wall" source and point sink 
on all sites, and the squares are from the "full wall" source with a point sink on even sites only. 
In this plot, the symbol size at each point is proportional to the confidence level of the fit, on a 
scale where the symbols in the legend correspond to 50% confidence. For this quark mass we used 
Dmin = 12 for our quoted mass. 
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FIG. 18. The nucleon to rho mass ratio in quenched and full QCD. The squares, diamond and 
octagons are the a ~ 0.13 fm matched lattice runs with zero, two and three flavors respectively. 
The fancy diamonds below the other points are preliminary quenched points at a ~ 0.09 fm, and 
the burst lying above the trend is a coarse lattice three flavor run at a ~ 0.2 fm. The octagon at 
the left is the physical value, and the octagon at the right is the trivial infinite quark mass value. 
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FIG. 19. The nucleon to rho mass ratio in quenched QCD with conventional and improved 
Kogut-Susskind quark action. The diamonds, squares, octagons and pluses are from the conven- 
tional action with lattice spacings of 0.16, 0.12, 0.07 and 0.043 fm respectively. The decorated 
pluses are the improved action results with a ~ 0.13 fm, and the fancy diamonds are preliminary 
improved action results at a ~ 0.09 fm. 
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FIG. 20. The decuplet baryon to vector meson mass ratio at a as 0.13 fm. Squares are the 
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E. Nonzero momentum 

As a check on the quality of the dynamics in our simulations, we calculated the energies 
of a few nonzero momentum mesons, namely the Goldstone pion with momenta ^ (0, 0, 1) 
and ^ (0, 1, 1) and the ^ z ® ^ z rho with momentum ^ (0, 0, 1). We compare these energies 
to the ideal dispersion relation by tabulating the "speed of light" in Table |X|. 

2 E(k)-E(6) 
° k 2 

We see that the dispersion relation is generally very good, but with noticeable deviations 
from one for the heavier mesons. The results are similar in the quenched and full QCD 



(7) 
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calculations, as can be seen by comparing lines with the same valence quark mass in Table [X|. 

VI. CONCLUSION 

We have used simulations with three flavors of dynamical quarks and a quenched sim- 
ulation on lattices with matched lattice spacing and physical size to isolate effects of the 
sea quarks on the hadron spectrum and on the static quark potential. This was done with 
a Symanzik improved gauge action and an improved Kogut-Susskind quark action to make 
the effects of the nonzero lattice spacing as small as practical. Effects of the sea quarks 
are clearly visible in the static potential. The mesonic mass ratio J is much closer to the 
experimental value when dynamical quarks are included. The clq meson couples strongly to 
two meson states, as expected when sea quarks are included. 

Several aspects need further study. While the quenched pion mass has the expected 
form, we have not understood the dependence of the three flavor pion mass on the quark 
masses. It would be very nice to be able to extract excited state masses, especially in the 
avoided level crossing of the do, and to extract decay rates to be compared to experiment. 
Although we have minimized the lattice artifacts by using an improved action, an empirical 
investigation of these effects is necessary. We are beginning a series of simulations with 
a ~ 0.09 fm to investigate these effects. 
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TABLES 



am Ujd 1 am s 


1 n / 2 

10/5 


Uq 


res. 


e 


lats. 


a/n 


quenched 


8.00 


0.8879 


na 


na 


408 


0.3762(8) 


0.02 / na 


7.20 


0.8755 


1 x 10~ 4 


0.013 


370 


0.3745(14) 


0.40 / 0.40 


7.35 


0.8822 


2 x 10" 5 


0.03 


332 


0.3766(10) 


0.20 / 0.20 


7.15 


0.8787 


5 x 10~ 5 


0.03 


341 


0.3707(10) 


0.10 / 0.10 


6.96 


0.8739 


5 x 10~ 5 


0.03 


339 


0.3730(14) 


0.05 / 0.05 


6.85 


0.8707 


1 x 10~ 4 


0.02 


425 


0.3742(15) 


0.04 / 0.05 


6.83 


0.8702 


5 x 10~ 5 


0.02 


351 


0.3765(14) 


0.03 / 0.05 


6.81 


0.8696 


5 x 10~ 5 


0.02 


564 


0.3775(12) 


0.02 / 0.05 


6.79 


0.8688 


1 x 10~ 4 


0.0133 


484 


0.3775(12) 


0.01 / 0.05 


6.76 


0.8677 


1 x 10~ 4 


0.00667 


407 


0.3852(14) 



TABLE I. Parameters of the a = 0.13 fm simulations 







0.01 


0.02 


0.03 


0.04 


0.05 


0.1 


0.2 


0.4 


H 


fit 


















7T 


L0 


18 


18 


18 


18 


18 


18 


18 


18 


K 


1,0 


14 


14 


14 


14 










ss 


1,0 


14 


14 


14 


14 










a 


2,1 


4 


4 


5 


5 


6 


6 


6 


6 


P 


1,1 


6 


7 


8 


9 


9 


10 


12 


14 


K* 


1,1 
1,1 


00 00 


8 
8 


9 
9 


9 
9 










a\ 


2,1 


4 


4 


4 


5 


5 


5 


6 


7 


bi 


2,1 


4 


4 


4 


5 


5 


5 


6 


7 


N 


1,1 


7 


8 


9 


10 


12 


14 


16 


18 


A 


1,1 


7 


7 


8 


9 












1,1 


7 


7 


8 


9 










A 


1,1 


3 


4 


5 


6 


7 


8 


9 


10 



TABLE II. Minimum distances used in propagator fits. With the exception of the A, these 
hadrons are obtained from the "corner" source. The top row is the light quark mass. The second 
column shows the type of fit used, where the two numbers are the number of simple exponentials 
included and the number of oscillating contributions included. For example, a fit of type "2,1" 
would include two particles with one parity and one particle with the opposite parity. Hadrons 
with nondegenerate valence quarks, such as the K, were computed only for m Ut d < rn s . 



33 



(mivalence 




amps 


range 




conf. 


n a / \ 
0.4 (IT) 


OO 


1.45664(12) 


18-31 


13/12 


n on 

0.39 


no/ \ 
0.2 (7TJ 


oo 


0.90106(16) 


1 O O 1 

18-31 


1 /I /l O 

14/12 


n 00 

0.32 


n i / \ 
U.l (IT) 


OO 


0.65693(19) 


1 O O 1 

18-31 


1 A 1 1 O 

14/12 


n 00 

0.28 


n n C / \ 

0.05 (7TJ 


oo 


n /i /^n /lo/ on\ 

0.40043(20) 


1 O O 1 

18-31 


15/12 


n c 

0.25 


n no / \ 

0.03 (7TJ 


oo 


0.35825(10) 


1 O O 1 

18-31 


01 / 1 
21/12 


n nc 

0.05 


0.02 (7r) 


oo 


n on A A r\ t ~\ t\ 

0.29440(17) 


18-31 


23/12 


0.026 


n m / 
0.01 (7TJ 


oo 


U.zllU4(lb ) 


18-31 


Oft 1 ~i 

zo/ 12 


n ni 
0.01 


U.U3/U.U5 (A J 


oo 


U.41zbl(ly J 


1 A Q1 

14-31 


Oft /1 ft 

zb/ 16 


n nft 
U.U6 


U.UZ/U.U5 (A J 


oo 


U.ooOOU^ZU ) 


1/1 Q1 


0^ /1 ft 
zo/ 10 


U.U I 


n ni /n / z<^\ 
U.Ul/U.Uo (A ) 


oo 


U.OOoOZ^Zo J 


1 /I Q1 
14-oi 


00 /ift 
Zo/ 10 


n 10 
U.lo 


n no /.^A 
0.02 [TTj 


n no 

0.02 


n onoco/ooA 
U.oUzoo^zzJ 


18-31 


n /i /1 
9.4/ Iz 


n ft7 

0.6 r 


n no /n nc / 7y"\ 

0.02/0.05 (A) 


n no 

U.Uz 


n onooo/ocN 


14-31 


00/ lb 


U.UU7b 


n n K / ^ tt\ 
0.05 (ssj 


n no 
U.Uz 


U.4/bzo(ZOJ 


1/1 1 

14-31 


z r/io 


(J. 114/ 


n /I / \ 
0.4 (IT) 


n /i 

0.4 


i /( /~*noo / i t\ 

1.46932(17) 


18-31 


c 0/10 

5.3/12 


n n r 

0.95 


n o /—"\ 
0.2 ("7TJ 


n o 

0.2 


n nTnon/oc\ 

0.97930(25) 


18-31 


00/10 

0.3/12 


0.76 


n i / \ 
0.1 (7TJ 


n 1 

0.1 


0. 68332(24) 


1 O O 1 

18-31 


1 a 1 1 

16/12 


0.17 


n n C / \ 

0.05 (7TJ 


n n c 

0.05 


n /iO/ioo/oi\ 

U.454zz(zlJ 


1 O O 1 

18-31 


O/"* / 1 O 

ZD/ 12 


n ni 1 

0.011 


n n /I /.~\ 
0.04 (7TJ 


n n^i /n nc; 

0.04/0.05 


(J. 43507(27; 


1 O O 1 

18-31 


1 A /I O 

14/ Iz 


n 1 

0.31 


n no / \ 

0.03 (ir) 


n no /n nr 

0.03/0.05 


0.37787(18) 


1 O O 1 

18-31 


1 r /1 

15/12 


n c 

0.25 


n no / \ 

0.02 {TTj 


n no / n nc 

0.02/0.05 


0.31125(10) 


1 o 01 

18-31 


16/ 12 


n 1 

U.zl 


n m / \ 
0.01 (IT) 


n n 1 /n nc 

0.01/0.05 


n o o i 1 o i o o \ 

0.22440(22) 


1 O O 1 

18-31 


1 /I / 1 O 

14/12 


n 0^7 

0.27 


n r\ a I n nc / ts~\ 

0.04/0.05 (A J 


n n a 1 n n c 

0.04/0.05 


n Af\\ ai /ot\ 

0.40141(27) 


1 A 1 

14-31 


1 O / 1 f 

18/16 


n 1 

0.31 


c\ no /n ncr / Ty\ 

0.03/0.05 (A) 


n no /n nc 

0.03/0.05 


n /ion oh n\ 

0.43013(19) 


1 A 1 

14-31 


26/16 


n n r 

0.052 


0.02/0.05 A~) 


0.02/0.05 


0.40984 21) 


14-31 


19/16 


0.28 


0.01/0.05 (K) 


0.01/0.05 


0.38334(29) 


14-31 


25/16 


0.072 


0.05 (ss) 


0.04/0.05 


0.48659(27) 


14-31 


18/16 


0.35 


0.05 (ss) 


0.03/0.05 


0.48796(18) 


14-31 


28/16 


0.035 


0.05 (ss) 


0.02/0.05 


0.49009(20) 


14-31 


23/16 


0.12 


0.05 (ss) 


0.01/0.05 


0.49443(25) 


14-31 


19/16 


0.26 



TABLE III. Pseudoscalar meson masses. Quenched results are first, followed by the single 
two-flavor run, followed by the three flavor runs. The first column is the valence quark mass(es), 
and the second column the sea quark mass or masses. The particle name in the first column is 
intended as a mnemonic. Here 'V indicates valence quark mass equal to the lighter dynamical 
quarks, or degenerate in the quenched case. "-FiT" indicates one valence quark equal to the light 
dynamical quarks and one at about m s , while "ss" indicates two valence quarks with mass about 
m s , in a flavor nonsinglet state. The remaining columns are the hadron mass, the time range for 
the chosen fit, x 2 and number of degrees of freedom for the fit, and the confidence level of the fit. 
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Q^lvalence 




amy 


range 


X l D 


conf. 


n a I ^\ 

0.4 (p) 


OO 


1.544b (2 J 


1 A 11 

14-31 


ID /1 /I 

18/ 14 


0.19 


n o / \ 

0.2 (p) 


OO 


i nnnn /r \ 

1.0900(5] 


1 O O 1 

12-31 


19/16 


n 07 

0.27 


n i / \ 

0.1 (p) 


OO 


0.8443(11) 


i n oo 

10-28 


O 1 / 1 E 

21/15 


n 1 e 

0.15 


n n e / „\ 

0.05 (p) 


OO 


n 71 po/oi \ 

0.7168(21) 


n 0*7 

9-27 


i n /1 e 

10/ 15 


n 

0.8 


U.U<3 (pj 


OO 


n /'•ceo/ o^A 

U.bboo(zbJ 


o o o 


n /i o 

b.y/lz 


n oc 

U.8b 


U.Uz [pj 


OO 


U.b4^z(oUj 


7 01 

f-zl 


7/11 

7/11 


n q 
U.o 


n m / ^ 

0.01 (p) 


OO 


n ^n'7 , n/'^n\ 
U.bU /U(bUj 


D-lO 


i n /7 

iu/ r 


n 1 n 
U.ly 


U.Uo/U.Uo (A ) 


OO 


n fim n ( Qn\ 
U.bylU^oU J 


O O/l 

y-Z4 


1 Q /l O 

lo/lz 


n qq 
U.00 


n no /n nE. ( z^*\ 
U.Uz/U.UD (A J 


OO 


n coon ^ Qn\ 
U.OoZU^oU J 


Q HI 

o-Zl 


in /in 
1U/ 1U 


n /i 
U.4Z 


n m /n nE / z^*\ 
U.U1/U.U5 (A J 


OO 


n c^on/' /i n\ 
U.bboU(4UJ 


Q Ol 

o-zl 


ii /in 
11/1U 


n a 
U.o4 


n no / ^^ 
U.Uz 


n no 

0.02 


U.bUUy(zo) 


7 OQ 


17/10 
1 7/l6 


n 
U.z 


n no /n nE / 7>^*\ 

U.U2/U.U5 (A J 


( 1 fill 

U.Uz 


n /^eoo/ooN 

0.6532(23) 


O 07 

8-27 


00 / 1 c 

22/ lb 


n 1 /i 

U.14 


n n e / 
U.U5 {(£>) 


n no 
U.Uz 


0. (003(14) 


o on 
8-29 


zo/lo 


n 01 
U.zl 


0.4 (p) 


0.4 


1 e /rio /o\ 

1.5602(2) 


1/1 1 

14-31 


1 / 1 A 

12/14 


0.64 


0.2 (p) 


n o 

0.2 


1 inci / 7\ 

1.1051(7) 


1 O O 1 

L2-6L 


1 A / 1 a 

14/ 16 


n E£? 

0.56 


n i / \ 

0.1 (p) 


n 1 

0.1 


0.8620(9) 


1 n o 1 

10-31 


1 a /1 

16/18 


n e 7 

0.57 


n n e / „\ 

0.05 [p) 


0.05 


n 71 E A /1 7^ 

0.7154(17) 


n o i 

9-31 


11 /in 

11/19 


n 01 

0.91 


0.04 (p) 


r\ r\ A /n n e 

0.04/0.05 


0.6853(17) 


n oo 

9-28 


11 /1 

11/16 


n 

0.83 


0.03 (p) 


n no /n n r* 

0.03/0.05 


n A c\c\ / 1 /i \ 

0.6490(14) 


8-27 


19/16 


0.26 


n no / ~\ 

0.02 (p) 


n no /n nr 

0.02/0.05 


n /?i 1 o/i nA 

0.5113(19) 


7 OO 

7-22 


IO/IO 

13/ 12 


n 

0.38 


n n 1 ( \ 

0.01 (p) 


n ni /n n e 

0.01/0.05 


n E7on/on\ 

0.5730(30) 


CIO 

6-18 


7A/9 


n En 

0.59 


n n /i /n nE / 7>^*\ 

0.04/0.05 (A J 


n n a /n n e 

0.04/0.05 


n 7n/in/ooN 

0.7040(22) 


n on 

9-29 


15/17 


n *2i 

0.61 


n no /n nE / 7>^*\ 

0.03/0.05 (A j 


n no /n n e 

0.03/0.05 


0.6845(18) 


n on 

9-29 


Ol / 1 7 

21/17 


n 00 

0.22 


0.02/0.05 (K*) 


0.02/0.05 


0.6631(19) 


8-25 


15/14 


0.38 


0.01/0.05 (K*) 


0.01/0.05 


0.6485(28) 


8-23 


4.7/12 


0.97 


0.05 ((p) 


0.04/0.05 


0.7198(19) 


9-29 


13/17 


0.77 


0.05 (<p) 


0.03/0.05 


0.7174(13) 


9-29 


19/17 


0.35 


0.05 (0) 


0.02/0.05 


0.7152(11) 


8-29 


21/18 


0.28 


0.05 ((f)) 


0.01/0.05 


0.7194(11) 


8-28 


11/17 


0.84 


TABLE IV. 


Vector meson 


masses. The format 


is the same 


as Table 
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Q^lvalence 


a^T-sea 


amsc 


range 


•.,2 / r> 
X /^ 


conf. 


0.4 (a ) 


OO 


i Ton / a \ 

1.739(4) 


/ ■ on 

6-20 


1 i /n 

11/9 


n on 

0.29 


no / ^ \ 

0.2 (a ) 


oo 


i one f£}\ 

1.29o[o) 


f on 
O-zU 


£ /n 
5/9 


n o o 


0.1 (aoj 


oo 


l.Ubo(loJ 


^ on 
o-zU 


1 o /n 
lz/9 


n o a 
U.z4 


U.Uo (clq) 


OO 


U.v (D{Zo ) 


D-ol 


oo /oo 
zz/zz 


C\ AQ 

U.4o 


U.Uo (clq) 


OO 


n n^nf a n\ 
U.yoU^4U J 


k on 
o-ZU 


o/in 

y/ iu 


n ^o 
U.oZ 


n no ( \n \ 
U.Uz l^oj 


OO 


U.o4o^Z f J 


a on 
4-ZU 


in/11 
1U/ 11 


n 


U.U1 (^oj 


OO 


n f!Qn/ /i n\ 
U.0oU(^4U ) 


a on 
4-ZU 


1 7 /I O 
1 // 1Z 


n 1 o 
U.1Z 


n no A 
U.Uz (do) 


n no 
U.Uz 


U.O /O^O J 


4-ol 


zo/ Z4 


n /1 1 
U.41 


0.4 (a ) 


n /l 

0.4 


1.750(4j 


£? on 

6-20 


o i /n 

21/9 


n m o 

0.013 


no/ \ 

0.2 (ao) 


r\ o 

0.2 


1.297(6) 


£? on 

6-20 


/i o /n 

4.8/9 


0.85 


0.1 (a ) 


0.1 


1.042(8 


6-20 


18/9 


0.031 


0.05 (a ) 


0.05 


0.829(6) 


5-31 


63/23 


1.3e-05 


0.04 (a ) 


0.04/0.05 


0.808(7) 


5-20 


10/10 


0.41 


0.03 (ao) 


0.03/0.05 


0.761(12) 


5-20 


16/10 


0.094 


0.02 (ao) 


0.02/0.05 


0.669(6) 


4-20 


31/11 


0.0011 


0.01 (a ) 


0.01/0.05 


0.532(8) 


4-20 


15/11 


0.2 



TABLE V. ++ (ao) meson masses. The format is the same as Table pll] . 



(MTlvalence 


a^sea 


ampv 


range 


X 2 /D 


conf. 


0.4 (ai) 


OO' 


1.816(7) 


7-20 


12/8 


0.17 


0.2 (ai) 


OO' 


1.370(7) 


6-20 


4.7/9 


0.86 


0.1 (ai) 


OO 


1.147(7) 


5-20 


3.6/10 


0.96 


0.05 (ai) 


OO 


1.037(13) 


5-20 


7.7/10 


0.65 


0.03 (ai) 


oo 


0.984(8) 


4-20 


15/11 


0.17 


0.02 (ai) 


oo 


0.950(12) 


4-18 


21/9 


0.013 


0.01 (ai) 


oo 


0.898(19) 


4-16 


21/7 


0.39 


0.02 (ai) 


0.02 


0.868(8) 


4-20 


11/11 


0.43 


0.4 (ai) 


0.4 


1.826(8) 


7-20 


12/8 


0.17 


0.2 (ai) 


0.2 


1.388(7) 


6-20 


12/9 


0.2 


0.1 ( ai ) 


0.1 


1.149(9) 


5-20 


7.3/10 


0.7 


0.05 (ai) 


0.05 


1.006(12) 


5-20 


28/10 


0.0018 


0.04 (ai) 


0.04/0.05 


0.964(10) 


5-20 


9.6/10 


0.47 


0.03 (ai) 


0.03/0.05 


0.925(6) 


4-20 


7.2/11 


0.78 


0.02 (ai) 


0.02/0.05 


0.874(6) 


4-20 


10/11 


0.53 


0.01 (ai) 


0.01/0.05 


8.131(11) 


4-15 


2.6/6 


0.86 




TABLE VI. 1 ++ (ai) 


meson masses. 


The format is the same 


as Table III. 
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Q^lvalence 




ampv 


range 


.,,2 / 7") 


conf. 


n a ( u \ 

0.4 (61) 


oo 


1 ooo/ n\ 

1.833(9) 


^7 on 

7-20 


13/8 


O 1 1 

0.11 


n o (u \ 

0.2 [pi) 


oo 




c on 
D-ZU 


o i /n 

8.1/9 


U.oo 


U.l (bi) 


oo 


1 1 CO/ 1 o\ 

l.loz(lo ) 


k on 
5-zU 


1 o /i n 
12/10 


n on 
U.zy 


U.UO \0\) 


OC' 


1 1 1 fl/ QC\\ 

l.llU^oU J 


on 


i o /i n 
lz/ 1U 


n q 
U.o 


U.Uo (^DiJ 


oo 


l.\j\J fylo) 


/i on 
4-ZU 


1C /I 1 

10/ 11 


n i ^ 
U.lO 


U.UZ (^Oij 


oo 


U.yoy^ZZ ) 


/i on 
4-ZU 


1 Q /I 1 

16/ 11 


n oc 
U.Zo 


U.U1 (^Dij 


oo 


n nnn/ c^n^i 
U.yyU^oU ) 


4- ID 


on/7 

o.y/ < 


n Oi^ 
U.Zo 


n no (u \ 
U.UZ (^ij 


n no 
U.Uz 


U.yUo(^14J 


/i on 
4-ZU 


11 / 1 1 
11/11 


n /i q 
U.4o 


0.4 (&i) 


0.4 


1.834(12] 


7-20 


1 o / o 

10/8 


n o / ' 

0.2b 


n o / l \ 

0.2 (bi) 


n o 

0.2 


1.398(10) 


n on 

6-20 


T Id 

7/9 


0.64 


0.1 (6i) 


0.1 


1.172(12) 


5-20 


7.9/10 


0.64 


0.05 (6i) 


0.05 


1.047(21) 


5-20 


2.5/10 


0.99 


0.04 (bi) 


0.04/0.05 


0.996(16) 


5-20 


11/10 


0.39 


0.03 (6i) 


0.03/0.05 


0.957(9) 


4-20 


9.5/11 


0.58 


0.02 (6i) 


0.02/0.05 


0.915(13) 


4-20 


13/11 


0.3 


0.01 (bi) 


0.01/0.05 


0.856(19) 


4-18 


8.7/9 


0.47 



TABLE VII. 1 + (b\) meson masses. The format is the same as Table 
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arriB 


range 


i n 
X IL> 


conf. 


n a /tvt\ 

U.4 (JNJ 


OO 


2.3933(15J 


18-31 


o 1 / 1 n 

21/10 


1 1 no 1 

U.U21 


n o /at\ 

U.2 (JNJ 


oo 


l.b847(23J 


lb-31 


i n / 1 o 

19/12 


n no /i 

U.U84 


U.l (JNj 


OO 


i o ^^n /rn\ 

1.277U(5UJ 


i A on 

14-30 


1 V 1 1 O 

15/13 


n o o 

U.33 


n n r /at\ 

U.U5 (JNJ 


oo 


l.U47U(7UJ 


1 O O 1 

12-21 


2.9/b 


n oo 

U.82 


n no /at\ 

U.U3 (JNJ 


oo 


U.94UU(5UJ 


n i n 

9-19 


8.1/7 


n oo 

U.32 


n no /at\ 

U.Uz (JNJ 


oo 


n ooon/ f 7n\ 

U.882U(7UJ 


O 1 T 

8-17 


e o la 

5.3/b 


n c 

U.o 


n m /at"\ 

U.U1 (JNJ 


oo 


n o 1 nn / nn\ 

U.81UU(9UJ 


T 1 A 

7-14 


£} O / A 

b.8/4 


n i c 

U.15 


n no /n nc / A \ 
U.Uo/U.Uo (AJ 


oo 


n o7Qn/ /i n^ 
U.y /oU(4U J 


q i n 

8-ly 


5. / /8 


n £ic 
U.b8 


n no /n / A , \ 
U.Uz/U.U5 (AJ 


oo 


n oqi n/ /i n^ 
U.yolU(4U J 


Tin 


t i /n 
/ .1/y 


n £io 
U.bo 


n m /n ne ( A , \ 
U.U1/U.U5 (AJ 


oo 


U.oyUU(bU J 


Tie 


11 / cc 
11/5 


n r\K 
U.Uo 


U.Uo/U.Uo (^J 


oo 


i nnnn/Qn 1 * 
!.UUyU(oU J 


q on 
o-ZU 


p. a /n 
o.4/y 


n 7 
U. 1 


n no /n nc / ,p zr > \ 
U.Uz/U.Uo (^J 


oo 


U.yoo / (zo J 


t on 
/-ZU 


c c /in 
0.0/ 1U 


n 

U. ( b 


n ni /n nc ( r zi\ 
U.Ul/U.Uo (^J 


oo 


U.yboU(4U J 


7 1Q 

f-lo 


n /i /c 
y.4/o 


n qi 
U.ol 


n no f~\T\ 
U.Uz (IN J 


n no 
U.U2 


r\ o a en/ /i n^ 
U.84oU(4U J 


q on 
8-zU 


i a, /n 

io/y 


n n^7/i 
U.Uf4 


n no /n ne / A \ 

U.Uz/U.Uo (AJ 


n no 

U.U2 


n m nn / on\ 

U.ylUU(oMJJ 


T 1 

7-21 


lb/11 


n 1 /i 

U.14 


n no /n nc^ 
U.Uz/U.Uo (^J 


n no 
U.U2 


n ott;o ( oq\ 
U.y /oz(zo J 


^7 OO 

l-ZZ 


1 //12 


n i k 
U.15 


U.4 (IN J 


n A 
U.4 


o /ioi q/i n 1 \ 
z.4zlo(lyJ 


1 O 91 

18-31 


/I /j /in 
4.0/ 1U 


n no 
U.92 


n o /at\ 

U.2 (JNJ 


n o 

U.2 


i 'ic\'~7 r - / o o \ 

1.7075(22 J 


lb-31 


17/12 


U.15 


n i /at\ 

U.l (JNJ 


n 1 

U.l 


1.311U(4UJ 


1 A O O 

14-28 


1)7 / l 1 

17/11 


n i o 

U.12 


U.U5 (JNJ 


U.Uo 


l.U57U(5UJ 


1 O O £ 

lz-zo 


e /i n 

5/1U 


n on 

U.89 


n n/i /at\ 

U.U4 (JNJ 


n n /i /n n cr 

U.U4/U.U5 


i nnon / on\ 

l.UU3U(3UJ 


1 n oo 

lU-zz 


i ^7 /n 

17/9 


f\ C\ A A 

U.U44 


n no /at\ 

U.U3 (JNJ 


n no /n n c 

U.U3/U.U5 


U.93UU(27J 


n oo 

9-22 


7.5/1U 


n no 

U.b8 


n no /at\ 

U.U2 (JNJ 


n no /n ncr 

U.U2/U.U5 


n o e /I n / on\ 

U.854U(3UJ 


o on 

8-2U 


c n /n 

b.9/9 


U.b5 


n m /at\ 

U.Ul (JNJ 


n n 1 /n nc 

U.Ul/U.Uo 


n '7'7C\C\{ (~T\\ 

U.779U(oUJ 


T 1 C 

7-lb 


1.1/0 


n no 

U.98 


n n /i /n n cr / A \ 

U.U4/U.U5 (AJ 


n n /i /n n cr 

U.U4/U.U5 


i no a n / on \ 

l.U24U(3UJ 


n oo 


n /i n 

9/1U 


n c o 

U.53 


n no /n ne / A \ 

U.U3/U.U5 (AJ 


n no /n n cr 

U.U3/U.U5 


U.9784(23J 


O 1 

8-21 


t /i n 

7/10 


n to 

U.72 


U U2/U U5 (A) 


U U2/U U5 


U 9312(27) 


7-20 


9.9/10 


0.45 


U.U1/U.U5 (A) 


U.U1/U.U5 


U.8850(5U) 


7-18 


4.2/8 


0.84 


U.U4/U.U5 (3) 


U.U4/U.U5 


1.0440(30) 


9-22 


8.3/10 


O.b 


U.U3/U.U5 (3) 


U.U3/U.U5 


1.0214(20) 


8-22 


8.3/11 


0.b8 


U.U2/U.U5 (3) 


U.U2/U.U5 


0.9989(20) 


7-21 


9.4/11 


U.58 


U.U1/U.U5 (3) 


U.U1/U.U5 


0.9798(27) 


7-2U 


6.5/10 


U.77 


TABLE VIII. Octet baryon 


masses. The format 


is the same 


as Table [TO]. 
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Q^lvalence 






range 




,,2 / n 

x /t> 


conf. 


n a ( A \ 

U.4 (A) 


OO 


O A O A / 1 \ 

2.424(1) 


10-25 




OO / 1 o 

39/12 


0.0001 


U.z 


oo 


1.737(2) 


9-25 




OO /l o 

30/13 


one 

0.0053 


0.1 (A) 


OO 


1.3M(b) 


Q 1 7 

8-1 / 




b.y/b 


n qq 
0.33 


u.uo 


oo 


l.lo4(lZ J 


l-Lo 




n Q7 /Q 

U.o / /o 


n n^ 
u.yo 


n HQ / A \ 

U.Uo (/a) 


oo 


i nno / 1 o\ 
l.Uyz(lz J 


C 1 1 

0-11 




o.o/o 


n qo 
U.oz 


u.uz 


OO 


l.UOo^lu J 


A 0. 

4-y 




U.oD/ z 


U.OD 


U.U1 \L\) 


OO 


1.UU<5(1 < J 


Q 7 




1 7/1 

1.7/1 


n in 

u.iy 


n no / a \ 
U.U2 (A) 


n no 

U.U2 


n nm /i o\ 

u.yyi(io ) 


/I i n 
4-1U 




1 T /O 
1. / /o 


n p. k 
O.bo 


U.4 (A) 


n A 

U.4 


2.450(b) 


1 O O CT 

18-25 




0.2/4: 


OIO 

0.19 


n o / A \ 

U.z (A) 


U.2 


1.753(2) 


o oo 

9-23 




12/11 


O OO 

0.33 


U 1 f A) 


U.I 


1 396(6) 


8-17 




11/6 


0.08 


U.U5 (A) 


U.U5 


1.155(8) 


7-15 




6.4/5 


0.27 


U.04 (A) 


U.U4/U.U5 


1.102(9) 


6-14 




8.9/5 


0.11 


U.03 (A) 


U.U3/U.U5 


1.059(6) 


5-13 




5.8/5 


0.32 


U.02 (A) 


U.U2/U.U5 


0.989(7) 


4-11 




10/4 


0.04 


U.01 (A) 


U.U1/U.U5 


0.921(12) 


3-9 




8.8/3 


0.032 


TABLE IX. Decuplet baryon masses. The format is the same as Table pTT[ 


dni'valence 




^(0,0,1) 




c*(o,m: 




c P (0,0,l) 


U.1U 


oo 


0.981(6) 




0.971(5) 




0.984(49) 


U.02 


oo 


1.001(10) 




0.973(21) 




0.900(82) 


U.U1 


oo 


0.996(14) 




0.981(59) 




1.021(122) 


U.U2 


U.U2 


0.997(11) 




0.986(16) 




1.013(62) 


U.4 


U.4 


0.966(19) 




na 




0.949(16) 


U.2 


U.2 


0.958(8) 




0.952(4) 




0.921(30) 


U.I 


U.I 


0.982(5) 




0.976(5) 




0.954(32) 


U.U5 


U.U5 


0.996(6) 




0.988(9) 




0.952(44) 


U.U4 


U.U4/U.U5 


1.000(8) 




0.972(9) 




0.926(51) 


U.U3 


U.U3/U.U5 


0.988(5) 




0.978(9) 




0.903(28) 


U.U2 


U.U2/U.U5 


1.001(8) 




0.993(16) 




0.982(60) 


U.U1 


U.U1/U.U5 


0.995(17) 




0.986(30) 




0.967(63) 



TABLE X. "Speed of light" for the pion and rho at a ~ 0.13 fm. The three columns are for 
the Goldstone pion with momenta ^7 (0, 0, 1) and ^7 (0, 1, 1) and the j z (g) j z rho with momentum 
f (0,0,1). 
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